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A combinatorial approach to the Euler Charac- 
teristic 

In the main body of this paper, we introduced the Euler characteristic (EC) and 
calculated its mean value per site for percolation clusters on two dimensional 
lattices. The approach presented in the main text can be generalized to higher 
dimensional cell complexes where the EC is given by the alternating sum of the 
number of cells of different dimension. 

Here, we want to present a complementary definition not restricted to cell 
complexes and better s uited to study c ontinuum percolation. Following the 
combinatorial approach ( |Hadwiger 1955 1, we start with the family K. d of closed, 
bounded, and convex subsets A C 
JC d by 



< dim(A) < d, and define the EC on 



X(A) 



1, 4e /C rf \0, 
0, A = e K d 



(1) 



Next, we construct a family TZ d of spatial patterns P by forming finite unions 
of convex bodies Ai € IC d , i.e. P — U"=i -^i, an d promote the EC to a valuation 
on the convex ring lZ d by requiring additivity (Hadwiger 1955) 



X(A U B) := X(A) + X(B) - X(A n B). 



(2) 



The value of X{A U B) is well-defined, since A, B £ IC d implies A n B e JC d . 
Iteration of Equation |2] leads to 



x({jA l ) = J2(-i) k ~ 1 E X(A ti nA l2 ...nA lk ), 



(3) 



k=l 



i\<i?...<ik 



for A\,... ,A n € JC d . Hadwiger has shown that the combinatorial EC agrees 
with the Euler-Poincare characteristic from algebraic topology, when the latter 
is restricted to the convex ring lZ d . It is now straightforward to turn these 
formal concepts into practical tools for the study of percolation clusters. 



Figure 1: Clusters on the kagome lattice. Left: To each white vertex, we attach 
a face of the dual lattice. White faces are connected if they share a dual edge 
or vertex, as indicated by the thick black lines. The white pattern displayed 
consists of two clusters Pi and P% with size(Pi) = 7 and size(P 2 ) = 1. The EC 
of the pattern is given by X(P 1 U P 2 ) = X{Pi) + X{P 2 ) =0 + 1. Right: The 
MEC of white clusters for occupation probability p. 

At first, we want to reproduce the relations presented in the main text. To 
this end, consider site percolation on a finite planar 2-d lattice A with periodic 
boundary conditions and let each lattice vertex be white with probability p and 
black with probability p = 1 — To construct spatial patterns corresponding 
to the white vertices, we attach to each white vertex a face of the of the dual 
lattice A*, as illustrated in Figure [l] for the kagome lattice. The fc-cells C£ of 
A*, i.e. its vertices v* , edges e* and faces /* for k = 0,1,2, are convex bodies 
as introduced previously and the EC of union of cells can be calculated using 
Equation [3] A white cluster is now defined as a maximal collection of white 
faces /* with nonempty pairwise intersections. Since a face is closed, meaning 
its boundary df* is part of the face (df* = /* H df*), two faces are connected 
if they share an edge e* or only a vertex v* . For example, the configuration 
of white cells on the kagome lattice shown in Figure [I] forms a pattern of two 
clusters Pi and P 2 . 

In percolation theory, the occurrence of a pattern P = U™^/* € 1Z 2 is a 
probabilistic event. Therefore, X(P) is now a random variable. Its mean value 
per site of the lattice A is defined by 

^ = ^ ^( U "i / *)^' (4) 

where N is the number of vertices of A. To evaluate the average of the EC, 
we use the iterated additivity relation from Equation [3] with the Aj replaced by 
/*, i — 1, . . . , ri^. The task is simplified by writing a closed polyhedral face /* 
as the union of its disjoint open components (open interior /* together with the 

1 Within the combinatorial approach, the connectivities of white vertices are more natural 
than the connectivities of black vertices; we therefore start out with a white pattern and 
recover the black patterns by complementation later. 
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open edges e* and the vertices v* from the boundary df*), and adopting the 
definition pefj 119811 iLikos et all 11995 1 



X(C* k ) = (-l) fe , fc = 0,l,2, 



(5) 



for the EC of open fc-dimensional lattice cells. With this representation of a 
white pattern, we obtain 



,L 2 ,L 1 ,L Q 



i—1 i—1 i—1 

= n* 2 — nl + Uq = no — n\ + n 2 . (6) 



Here n* k , k 



0, 1, 2 denotes the number of white faces, edges and vertices of the 
pattern on A* and the last equality follows from duality. The average of the EC 
depends on the specific geometry of the underlying lattice. As an example, let 



us consider the eleven Archimedian lattices (comp. main text, (Griinbaum and 
Shephard, 1986[ |). In an Archimedian lattice with vertex type (ni, . . . ,n z ) and 
periodic boundary conditions, the fractional numbers of bonds and polygonal 
faces with m edges per site of the lattice A are given by | and ^- , respectively. 
With this structural information, we find the averages 



1 

N { 
1 



N 



(I-P 2 ), 



(7) 



These expressions are obtained by noting that a face /* is white with probability 
p = 1—p; an edge e* is white with probability 1 — p 2 , since it belongs to a white 
pattern if at least one of the faces sharing e* is white; likewise, a vertex v* is 
white, if at least one of the faces sharing v* is white, and this happens with 
probability 1 —p ni . Thus from Eq. [5] [6] and [7] we arrive at the result 



m = P - [i - a - pf] ph 



(8) 



for the MEC of percolation patterns on Archimedian lattices. For example, 
Equation [8] yields for the kagome lattice (3, 6, 3, 6) 

XkagiP) = p(l -p)(l-ip+ yp 2 - ^ 3 + ^ , (9) 

and for the triangular lattice 



Xtri{p) =P(1 -p) (1 - 2p) . 



(10) 
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A given white pattern P = U™^/* determines a unique black pattern P, as the 
set-complement, P com p, of P with respect to the lattice A, and we have 

X(P(JP)=X(P) + X(P)=0. (11) 

The first equality holds because PP\ P = by definition of the complement, the 
second equality follows since P U P is topologically a two dimensional torus. 

Obviously, the complementary pattern P may be viewed as an aggregate of 
percolation clusters of faces of A* chosen with probability p = 1 — p and colored 
black. The MEC of black patterns is defined by 

X (p)= lim hx{P)) p . (12) 

TV— >oo iv 

Consequently, Equation 1 1 1 1 implies the symmetry-type relation 

X(1-P)+X(P)=0. (13) 

Thus, the MEC of black patterns is obtained from that of white patterns via 
simple substitution. After substitution, we obtain exactly the same polynomials 
for xip) as m the main text, for example: 

Xkag(p) = P(l-p) [1-P-\{P 2 +P 3 +P^) 



3 

Xtriip) = p(l - p)(l - 2p) = Xtri(p) (14) 

We note, that the trianglular lattice is the only one among the Archimedian 
lattice with equal "black" and "white" MECs. 
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Summary of two-dimensional lattices 



In this section, we collect all formulae for the mean Euler characteristic of the 
various lattices discussed in the main text and present tables containing the 
numerical values of po, p* , and p c . We also discuss some irregular lattices 
and the dual Archimedean lattices, which are not included in the main paper. 
Maple™ sheets used to calculate the MECs as well as numerical values of po 
and p* are available from the authors on request. 



Archimedean lattices — site percolation 

Each vertex of an Archimedean lattice is surrounded by z regular polygons with 
ni, . . . , n z edges. In terms of the n^, the MEC of an Archimedian lattice is given 
by 



/ Z 1 7li— 3 \ 

x(p)=p(i-p) i-pE- EH 

\ i=l n% fj.=0 J 



(15) 



The vertex configurations, and a comparison of the zero crossing, our estimator 
of Pc and the actual threshold are given in Table [I] 



m, ■ .-,n z 


Po 


P 


pSite 


3,12,12 


0.8395 


0.7869 


0.807900 . . . 


4,6,12 


0.7833 


0.7373 


0.747806(4) 


4,8,8 


0.7689 


0.7269 


0.729724(3) 


6,6,6 


0.7413 


0.7043 


0.697043(3) 


3, 6, 3, 6 


0.6756 


0.6462 


0.6527036 . . . 


3,4,6,4 


0.6468 


0.6224 


0.621819(3) 


4,4,4,4 


0.6180 


0.5987 


0.5927460(5) 


3, 3, 3, 3, 6 


0.5913 


0.5752 


0.579498(3) 


3,3,4,3,4 


0.5616 


0.5511 


0.550806(3) 


3,3,3,4,4 


0.5616 


0.5511 


0.550213(3) 


3, 3^ 3j 3 ; 3, 3 


0.5 


0.5 


0.5 



Table 1: Archimedean lattices - site percolation. Lattices are ordered with 
decreasing p c , the percolation thresholds are from Suding and Ziff (19991 and 



references therein. These values are plotted in Figure 3 of the main paper. 



Archimedean lattices — bond percolation 

The MEC of the covering lattice of an Archimedean lattice with vertex type 
tlx , . . . , n z is given by 

x(p) = -p+ 2 (i-(i-p) z ) + Y,— p rH - ( 16 ) 
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ni, ... ,n z 


Po 


P 


rj bo7id 
I c 


3,12,12 


0.7580 


0.7098 


0.7404219(8) 


4,6,12 


0.7054 


0.6685 


0.6937338(7) 


4,8,8 


0.6964 


0.6623 


0.6768023(6) 


6,6,6 


0.6756 


0.6462 


0.6527036 . . . 


3, 6, 3, 6 


0.5277 


0.5227 


0.5244053(3) 


3,4,6,4 


0.5134 


0.5111 


0.5248325(5) 


4,4,4,4 


0.5 


0.5 


0.5 


3, 3, 3, 3, 6 


0.4069 


0.4233 


0.4343062(5) 


3,3,4,3,4 


0.3992 


0.4166 


0.4141374(5) 


3,3,3,4,4 


0.3992 


0.4166 


0.4196419(4) 


3} 3} 3; 3} 3j 3 


0.3244 


0.3538 


0.3472963 . . . 



Table 2: Archimedean lattices - bond percolation. Lattices are in the same order 
as in Table [T] numerical estimates of percolation thresholds are from |Parviainen| 
(20071. This data is plotted in Figure 3 of the main paper. 



Dual Archimedean lattices — site percolation 



The MEC of the dual lattice of an Archimedean lattice depends only on the 
coordination number of the Archimedean lattice. 



xip) = p0--p) 1 




(17) 



A graph illustrating the relationship between p c and the zero crossing of the 
MEC and the performance of our estimator p* is shown in Figure [2j The dual 
Archimedean lattices arc also known as Laves lattices. 



dual(ni, . . . , n z ) 


Po 


* 

P 




dual(3, 12, 12) 


0.5 


0.5 


0.5 


dual(4, 6, 12) 


0.5 


0.5 


0.5 


dual(4,8,8) 


0.5 


0.5 


0.5 


dual(6, 6, 6) 


0.5 


0.5 


0.5 


dual(3, 6, 3, 6) 


0.6180 


0.5987 


0.5848(2) 


dual(3,4,6,4) 


0.6180 


0.5987 


0.5824(2) 


dual(4,4,4,4) 


0.6180 


0.5987 


0.5927460(5) 


dual(3, 3, 3, 3, 6) 


0.6914 


0.6612 


0.6396(2) 


dual(3,3,4,3,4) 


0.6914 


0.6612 


0.6500(2) 


dual(3,3,3,4,4) 


0.6914 


0.6612 


0.6476(2) 


dual(3, 3, 3, 3, 3, 3) 


0.7413 


0.7043 


0.697043(3) 



Table 3: Dual Archimedean lattices - site percolation. Lattices are in the same 
order as in Tab le [Tj numerical estimates of percolation thresholds are from 
20031; see Figure [2] for a plot of this data. 



van der Marck 
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2-uniform lattices — site percolation 



While Archimedean lattice have only one vertex type, the vertices of 2-uniform 
lattices are of two distinct types. A vertex is again characterized by the polygons 
surrounding it. Each vertex type i makes up a fraction Sj of the total vertices 



and a 2-uniform lattice is typically denoted by si(n\, . . . , n\ ) + S2(nf , . . 
The MEC can be expressed in terms of the vertex configuration: 



— 3 



X(p) =p(l-p) 



v=l,2 i=l 1 fj,=Q 



J- 



(18) 



We determined the percolation thresholds of all 2-uniform lattices using the 
algorithm by Newman and Ziff (Newman and Ziff 2001). Simulation were 



performed on lattices with linear dimension L — 64, 128, 256, 512 and 1024. The 
probabilities at which a spanning appears on a finite lattice were extrapolated 
to infinite lattice size using the well known finite size scaling relations. 



(n\, . . .,nlj + s 2 (n(, . 



3,4,3,12) + £(3,12*) 
3, 4, 6, 4) + |(4, 6, 12) 
3 6 ) + f (3 2 ,4, 12) 
3 2 ,6 2 ) + i(3,6,3,6) 
3 6 ) + §(3 2 ,6) 
3,4 2 ,6) + i(3,6,3,6) 
3,4 2 ,6) + ±(3,6,3,6) 
3,4 2 ,6) + |(3,4,6,4) 
3 4 ,6) + ±(3 2 ,6 2 ) 
3 3 ,4 2 ) + ±(3,4,6,4) 
3 2 ,4,3,4) + i(3,4 
3 3 ,4 2 ) + i(4 4 ) 
3 3 ,4 2 ) + i(4 4 ) 
3 6 ) + |(3 4 ,6) 
3 3 ,4 2 ) + i(3 2 , 4,3,4) 
3 3 ,4 2 ) + |(3 2 ,4,3,4) 
3 6 ) + f(3 2 , 4,3,4) 
3 6 ) + i(3 4 ,6) 
3 6 ) + |(3 3 ,4 2 ) 
3 6 ) + |(3 3 ,4 2 ) 



Po 



P 



6,4) 



0.7909 
0.7424 
0.6918 
0.6756 
0.6532 
0.6526 
0.6526 
0.6468 
0.6354 
0.6053 
0.6053 
0.5907 
0.5811 
0.5684 
0.5616 
0.5616 
0.5530 
0.5453 
0.5414 
0.5311 



0.7403 
0.7016 
0.6555 
0.6462 
0.6270 
0.6271 
0.6271 
0.6224 
0.6119 
0.5874 
0.5874 
0.5755 
0.5675 
0.5563 
0.5511 
0.5511 
0.5439 
0.5374 
0.5343 
0.5258 



0.7680(2 
0.7157(2 
0.6808(2 
0.6499(2 
0.6329(2 
0.6286(2 
0.6279(2 
0.6221(2 
0.6171(2 
0.5885(2 
0.5883(2 
0.5720(2 
0.5648(2 
0.5607(2 
0.5505(2 
0.5504(2 
0.5440(2 
0.5407(2 
0.5342(2 
0.5258(2 



Table 4: 2-uniform lattices site percolation. Lattices are ordered according to 
decreasing percolation threshold. Compare with Figure 3 of the main paper. 



Irregular lattices 

For lattices that are not included in any of the above classes classes, we present 
the polynomials x(p) here. Table [5] and Figure [3] summarize our results for 
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site percolation: Laves lattices 




5 10 

Figure 2: Dual Archimedean lattices - site percolation. The solution of Equation 
18 of the main text, p* , is a good estimate of p s ^ te for all dual Archimedian 
lattices, also known as Laves lattices. For numerical values of p a , p* and p c , see 
Table El 



irregular lattices. 

In addition to the Archimedian lattices, 



a number of other lattices have been 



studied by Suding and Ziff (1999), including the bowtie lattice with vertex con- 
(3 4 , 4 2 ) + | (3 2 , 4 2 ) the dual bowtie lattice with vertex configuration 



figuration | 

^ (4 2 ,6 2 ) + | (4,6 2 ). Given these vertex configurations, x(p) can be calculated 
using the same formula as for 2-uniform lattices. 



Voronoi tesselation and its dual. Consider points randomly and indepen- 
dently distributed on the plane and construct faces around these points using 
the Wigncr-Scitz construction. The resulting Voronoi tesselation of the plane 
consists of faces each of which contains exactly one of the randomly distributed 
points. The faces have six sides on average and the coordination number of 
a vertex is three with probability one. The dual of the Voronoi tesselation is 
a completely triangulated lattice with the randomly distributed points as ver- 
tices. The MEC of the dual tesselation is therefore the same as that of the 
triangular lattice. To calculate the MEC of the Voronoi lattice, we need to 
know the distribution Vi of the faces with different number of sides. This is 



known approximately (D. Stoyan 1996) and the MEC is found to be 



Xvorip) = p - 1 + -(1 - p 2 ) - ^X^ 1 



(19) 



i>3 



We are not aware of any estimation of the site percolation threshold for the 
Voronoi tesselation. 
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0.8 



0.7 



0.6 



0.5 



Voronoi 




Delauney 







1 
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Figure 3: Irregular lattices - site percolation. The relationship between po, 
p* and p s ^ te also holds for irregular, random and quasi-periodic lattices. For 
numerical values, see Table [5] 



lattice 


si(nj, . . . ,<J + s 2 (n(, . . . ,n£J 


Po 


P* 


l J C 


bowtie 


A (3, 3, 3, 3, 4, 4) + A (3, 3, 4, 4) 


0.5616 


0.5511 


0.5474(8) 


bowtie dual 


A (4, 4, 6, 6) + |(4, 6, 6) 


0.7048 


0.6722 


0.6649(3) 


dice 


A (4, 4, 4, 4, 4, 4) + 2(4, 4, 4) 


0.6180 


0.5987 


0.5848(2) 


Penrose 




0.6180 


0.5987 


0.5837(2) 


Penrose dual 




0.6602 


0.6334 


0.6381(3) 


Voronoi 




0.7548 


0.7151 




Voronoi dual 




0.5 


0.5 


0.5 



Table 5: Irregular lattices 
from 



(van der Marck 
Figure[3j 



site percolation. Percolation thresholds are taken 
1997 Yonezawa et al. 1988). This data is plotted in 



Penrose tiling and its dual. The rhombic Penrose tiling is a quasi periodic 
lattice consisting of 4-sided polygons only. The vertices have an average coordi- 
nation number of 4, varying between 3 and 7. The distribution of the different 



vertex types is known (Lu and Birman 1987), but since the tessellation consists 
of 4-sided polygons only, the MEC is the same as that of the square lattice. 
The dual Penrose tiling, however, has faces with 3 to 7 sides. The coordination 
number is uniformly z — 4, and from the distribution of faces we find for the 
MEC 

Xip) = -(1 - p) + 2 (1 - p 2 ) - Vi (1 - P l ) , (20) 

i>3 



where z/j is the probability that a given face has i sides (Lu and Birman 19871. 
The percolation thresholds, po and p* for these lattices are given in Table]5j 
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Three dimensional cubic lattices 



The generalization of Equation [5] to three dimensions states that the MEC of 
the white pattern is given by the alternating sum of the mean number of cells of 
different dimensions per vertex. For cubic lattices, natural space filling cells are 
given by the Wigner-Seitz cells. While white vertices are connected whenever 
their cells share a face, an edge or a vertex, black vertices are connected only 
when their cells share a face. For the simple cubic lattice, the MEC is given by 

XsM=P~'3p 2 + 3 P 4 -p 8 . (21) 

The MEC of the face centered cubic (fee) lattice is given by 

X fcc ( P ) = P (l-p)(p 4 +p 3 + 3p 2 -5 P +l). (22) 

A vertex of the fcc-lattice has 12 nearest neighbors, each of which corresponds 
to one face of the WSC. 

A vertex of the bcc-lattice is connected to eight neighbors. The Wigner-Seitz 
cell, however, has 14 faces, six of which make contact to next nearest neighbors 
on the lattice. The MEC obtained when using the connectivities of the Wigner- 
Seitz cells corresponds to a lattice where the six next nearest neighbors are 
adjacent. 

Xu-uip) = P(l - P)(6p 2 - 6p+ 1) (23) 

The MEC of the bcc-lattice for black and white vertices are identical, reflecting 
the fact that at any 1-cell (edge) or a 0-cell (vertex) of the WSC no cells meet 
that do not share a face. 
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